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INTRODUCTION 


Clear air turbulence (CAT) encountered by aircraft in the 
troposphere and ietstream regions extending above the tropopause 
has been empirically related to numerous meteorological features 
and synoptic patterns. Several studies (refs. 1 to 9) indicate 
that the principal physical features which are critical to the 
occurrence of CAT include vertical wind shear and atmospheric 
stability. As higher altitude aircraft have been developed 
and have encountered turbulence in the stratosphere (refs. 10 
to 14), the examination of meteorological conditions has also 
come to include vertical wind shear and static atmospheric 
stability (refs. 15 to 20). The data available for the study 
of high altitude turbulence (HAT) inherently include localized 
microscale and mesoscale wind shear and lapse rate features 
for which no adequate climatology has been established. 

Moreover, a suitable climatology cannot be based on meteoro- 
logical data acquired to date, since HAT projects have typically 
been limited to a few case studies and the meteorological 
elements of interest have varied with the concerns of each 
investigator. In addition, data for some case studies resulted 
from deliberate attempts to encounter HAT, Thus, without an 
adequate climatology, any interpretation of the statistical 
relationship between meteorological features and HAT must be 
guarded. For example, one must inquire whether the threshold 
levels for individual variables more accurately specify turbu- 
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lence, nonturbulence, or both. Further, it is not clear 
whether the association between individual meteorological 
variables and HAT is due to physical cause-effect relationships 
or simply to coincidental seasonal patterns. This problem is 
of greater significance to the practitioner who wishes to 
apply the results of individual research studies to operational 
forecasting procedures than it is to the scientist who wishes 
to identify the physical processes that generated a particular 
turbulence event. 

The study discussed herein was initiated to obtain prelim- 
inary indications of the climatology of selected variables 
that were expected to be associated with the occurrence or 
nonoccurrence of high altitude turbulence. Both theoretical 
considerations and empirical evidence show that certain combin- 
ations of vertical wind shear and temperature lapse rate are 
associated with turbulence. The classical concept of Richardson 
number (Ri), the ratio of thermal stability to the square of 
the vertical wind shear, indicates that turbulence will prevail 
or propagate when Ri is less than 1; that is, when the kinetic 
energy available from wind shear exceeds the energy needed to 
overcome the static stability or buoyancy forces in the atmos- 
phere. In addition, when dynamic conditions favor Kelvin- 
Helmholtz instability, gravity wave motion can amplify and 
subsequently break into turbulence. The critical value of Ri 
for this process has been shown to be less than 0.25 
(refs. 2 and 21). This process has been observed to result in 
strong turbulence when wind shear increases to large magnitudes 
in the presence of high static stability (refs. 3, 7, and 15). 
The primary variables chosen for the present study are the 
strongest positive lapse rates, negative lapse rates, vertical 
wind shear, and wind shear windspeed product (vertical gradient 
of kinetic energy) that occur within separate layers of the 
upper troposphere and lower stratosphere. Wind shear in the 
troposphere and windspeed at selected mandatory levels are 
also included. This paper reports some of the initial frequency 
distributions and data problems. Examples of the application 
of the climatological data are discussed, and recommendations 
for further study are made. 


SYMBOLS AND ABBREVIATIONS 


CAT clear air turbulence 

HAT 
HI CAT 
n 


2 


high altitude turbulence 

high altitude clear air turbulence 

number of observations in a statistical sample 


P( ) 






Ri 



V 


probability of an event ( ); the conditional proba- 
bility of event y based on knowing that event x 
does occur is designated P(y/x) 

Richardson number, the ratio of the energy needed 
to overcome atmospheric buoyancy to the kinetic 
energy available from vertical wind shear 

variance of an observed variable 

turbulence intensity categories used in reference 17, 
nominally characterized as negligible, light, and 
moderate, respectively 

wind velocity in the horizontal plane 


AZ 


vertical wind shear, the rate of change of the hori- 
zontal wind velocity with altitude 


X 


dummy variable 


Z vertical Cartesian coordinate, also designation for 

hours (and minutes) in Greenwich mean time 

z normal distribution standardized variable equal to 


rt measure of uncertainty; the confidence level equals 

1 - It 

Y lapse rate for a layer of low static stability 

( late of ambient temperature decrease with altitude), 
"C/kro 

-Y temperature inversion rate for a layer of high static 

stability, "Ckm 

A incremental operator 

p mean value 

,1^ true vaiiance of the population tiom which the 

variable is sampled 

chi-square distribution variable 
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Subscripts: 

fit 


flight level 
low altitude 
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low 


max 

maximum 


100 

100 millibar 

level 

300 

300 millibar 

level 

Seasons: 



w 

winter (Dec. 

to Feb . ) 

SP 

spring (Mar. 

to May) 

su 

summer (June 

to Aug . ) 

F 

fall (Sept. 

to Nov. ) 

A 

annual (all 

seasons combined) 

Turbulence 

events : 


NF 

not forecast 


NR 

not encountered (not recorded) 

TF 

forecast 


TR 

encountered 



DATA AND PROCEDURES 


Since this was a pilot study, the data sample, class 
intervals, and altitude layers of interest were selected 
somewhat arbitrarily, with emphasis on phenomena above the 
troposphere. Consideration was given to the quality of the 
initial upper air observations, the results of turbulence 
studies, and the suitability of the data for automatic pro- 
cessing. For the period from July 1, 1970 to June 30, 1971, 
the 002 and 122 radiosonde observations were examined for 20 
upper-air stations in the western United states. This period 
was selected because the Wi.nd and thermal upper air data taken 
after July 1, 1970 were consolidated into new and much improved 
archival formats at the National Climatic Center. In addition, 
more rigorous quality control procedures were implemented at 
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' this time. The 20 stations selected were San Diego, Oakland, 

and Vandenberg Air Force Base, Calif.; Yucca Flat, Ely, and 
Winnemucca, Nev.; Denver and Grand Junction, Colo.; 

■ Albuquerque, N. Mex.; Tucson and Winslow, Ariz.; Lander, Wyo.; 

Salt Lake City, Utah; Boise, Idaho; Great Falls and Glasgow, 

^ Mont.; Spokane and Quillayute, Wash.; and Medford and Salem, 

( oreg. The height, temperature, and both zonal and meridional 

wind components were recorded for each mandatory and significant 
, pressure level from 700 millibars to the termination of the 

I sounding. 

Table 1 lists the class intervals and restrictions for 
the variables examined. The maximum values of each variable 
I within the indicated layer were summarized seasonally and 

) annually after the data for both observation times for all 

I stations were combined. Univariate percent frequencies and 

I cumulative percent frequencies were computed for all variables 

) for all layers. Bivariate percent frequencies were also 

computed for selected elements and layers. 

The computational procedures involved the partitioning of 
the layers given in table 1 into sublayers defined by mandatory 
and significant levels in the radiosonde observations. Meteor- 
' ological variables were computed for overlapping sublayers 

I after insuring that the restrictions and threshold values in 

table 1 were met. The maximum value of the variable in the 
given layer was extracted and retained for later analysis. It 
is important to note that for this maximum value, the values 
of other variables that occurred simultaneously were not 
^®f^tned. For example, if the maximum wind shear was observed 
in the sublayer between 135 and 150 millibars, the corresponding 
values in this sublayer of lapse rate and the wind shear 
' windspeed product were not retained in the data base unless 

they also met their respective threshold criteria. This 
deficiency presented serious problems in the subsequent determi- 
nation of multivariate relationships. 

The values of the tabulated univariate frequency data for 
the parameters listed in table 1 were plotted against cumulative 
percent frequencies on normal probability versus logarithmic 
scales. Both annual and seasonal data for each layer were 
gtaphed for most of the variables; the graphs are pi'esented in 
figures 1 through 67 , which are discussed in the next section. 

The lines of fit are manual fairings that are intended to show 
gross features of the distributions and not lines of best fit 
based on least squares or other theory. The threshold values 
given in table 2 are represented by thin vertical lines on the 
graphs. A logaritlimic scale was chosen because the lower 
magnitude class intervals contained most of the observations, 
thus skewing the frequency distributions. Using this scale, a 
straight line fit would indicate a log-normal probability 
density function. 
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The authors recognize that probability plotting and 
manual fairing is subjective in that the validity of any 
assumed statistical model is based on visual examination. The 
advantages of the technigue are that the method is simple, the 
data are pictorially represented, the reasonableness of an 
assumed model is easily evaluated, and the parameters and 
percentiles of the distribution can be estimated. The main 
disadvantages are the lack of both objectivity and of a proba- 
bilistic framework that could be provided by statistical 
calculations and tests for the distribution functions. 


RESULTS 


Univariate Distributions 

Results are presented for each variable in the sequence 
of annual and then seasonal frequency distributions. Data for 
wind shear and lapse rate are given first for combined layers 
and then for the individual layers in the lower stratosphere. 

The distributions of windspeed with altitude are shown in 
figures 1 to 5 for the data sample used in this study. The 
300 millibar level shows the highest windspeeds and is consistent 
with the midlatitude jetstream that occurs just below the 
tropopause. Speeds general Iv decrease up to about 50 millibars 
and then increase slightly. The increase in windspeeds above 
50 millibars (70 mb in summer) is associated with the change 
in circulation patterns between the troposphere and strat- 
osphere. Tropospheric midlatitude winds are predominantly 
westerly throughout the year, while lower stratospheric winds 
are predominantly westerly in winter and easterly in summer. 

This stratospheric wind reversal is strongly reflected in the 
seasonal windspeed distributions at 50 and 30 millibais. At 
each level through 70 millibars (figs. 6 to 9) the windspeeds 
are highest in winter, lowest in summer, and slightly higher 
in spring than in fall. At 50 millibars (fig. 10) the winds 
are strongest in winter and about the same through the othei 
rhree seasons. At 30 millibars the influence of the shift 
from westerlies to easterlies causes the summer windspeeds to 
exceed those in winter to the 75 percentile frequency of 
occurrence (fig. 11). 

On either an annual or seasonal basis, the distributions 
for maximum wind shear within the lower tropospheie, 700 to 
350 millibars, and the lower stratosphere, 175 to 20 millibars, 
are similar (figs. 12 to 16). Stronger values of the maximum 
wind shear, 0.014 to 0.03 per second, occur relatively more 
often in the upper troposphere, 350 to 175 millibars. The 
distributions for these three broad atmospheric layers meige 
and change shape between shear values of 0.03 and 0.05 per second. 
This change in the distribution near the 99.9 percentile 
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suggests the presence of Measurement errors in the data. 
Seasonally, strong wind shears are most frequent during winter 
and least frequent during summer in each of these three layers 
(figs. 17 to 20). 

The frequency of strong wind shears in the individual 
layers above 175 millibars shows a consistent decrease with 
altitude from the 175 to 125 millibar layer to 40 to 20 millibar 
layer during each of the seasons as well as on an annual basis 
(figs. 21 to 25). As shown in figures 26 to 30, each of these 
layers also contains strong wind shears most frequently in 
winter and least frequently in summer. This seasonal pattern 
is most pronounced in the upper layers, in spite of the tendency 
for the windspeeds at these levels to be lowest in the spring 
rather than in the summer. The relative absence of strong 
wind shears in summer even though windspeeds are not at a 
minimum suggests a decreased occurrence of dynamic perturbations 
at 50 millibars and 30 millibars in summer. 

Distributions for the product of wind shear and windspeed 
(figs. 31 to 35) show decreasing magnitudes as altitude increases. 
^ typical seasonal behavior, with small magnitudes most prevalent 
in summer and large magnitudes most prevalent in winter, is 
observed in all layers except for the 40 to 20 millibar layer 
(figs. 36 to 41). In this layer, large magnitudes occur least 
frequently in the spring, when windspeeds are lowest, and in 
summer, when wind shears are minimal. 

The annual maximum positive lapse rate observed within a 
layer decreases in magnitude with increasing layer altitude, 
as shown in figure 42. This decrease with altitude is also 
characteristic of the individual seasons (figs. 43 to 46) 
except for the 175 to 125 millibar and 125 to 80 millibar 
layers. In winter and spring the magnitudes of the lapse 
rates do not decrease noticeably with altitude from the 175 to 
125 millibar layer to the 125 to 80 millibar layer but remain 
approximately the same. During these seasons the layers above 
175 millibars are usually in the lower stratosphere, where 
conditions are controlled by the intense midlatitude cyclonic 
activity in the troposphere and by the energy and momentum 
transfer through the tropopause. In general the seasonal 
patterns (figs. 47 to 54) show higher positive lapse rates in 
winter in the stratosphere and in summer in the troposphere. 

These patterns are a reflection of an increase of dynamic 
perturbations in the winter stratosphere and increased heating 
of the earth and convective activity in the troposphere during 
the summer. 

Above 125 millibars, strong inversions occur less frequently 
as altitude increases (figs. 55 to 59). In the lower layers 
strong inversions are more prevalent between 350 and 175 millibars 
than between 700 and 350 millibars. The 175 to 125 millibar 
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layer is similar to the 350 to 175 millibar curves in the 
summer and fall and to the 125 to 80 millibar curves in winter 
and spring. The correlation of the 175 to 125 millibar layer 
with tropospheric conditions in the summer or with stratospheric 
conditions in the winter is directly dependent upon the seasonal 
shift in the height of the tropopause. 


The seasonal patterns within layers (figs. 60 to 67) are 
variable. Between 700 and 350 millibars and between SO and 
40 millibars, strong inversions occur least often in summer 
and with similar frequency in the other three seasons. Between 
350 and 125 millibars they occur more often in winter and 
spring and less often in summer and fall. Virtually no seasonal 
pattern exists in the 125 to 80 millibar layer. In the 40 to 
20 millibar layer inversions are strongest in winter, weakest 
in summer, and about the same magnitude in spring and fall. 

The maximum inversions in the lower stratosphere as a whole 
(tropopause to 20 mb) are strongest in winter, weakest in 
summer, and stronger in spring than fall. 


Bivariate Distributions 

The joint occurrence of two variables with values of 
®*^f^icient magnitude for both to affect the development of 
turbulence is of particular significance to operational fore- 
casting procedures. If they frequently occur together, for 
example, and are highly correlated, then only one need be 
included in an objective scheme. On the other hand, if both 
variables affect the development of turbulence but are not 
correlated with each other, both should be included. 

Therefore, bivariate contingency tables were prepared to 
examine the dependence of positive and negative lapse rates 
and high altitude wind shears on the other variables, as well 
as on each other. No attempt was made to arrive at distribution 
assumptions. The bivariate tables were reduced to four-cell 
contingency tables by arbitrarily assigning threshold values 
to each variable above which noticeable high altitude turbulence 
might be expected. These values are shown in table 2. The 
technique estimates those variables that could be important in 
an objective forecasting scheme such ar a multivariate linear 
regression model . 

A difficulty in the use of the present data set was 
experienced when an objective chi-square test was applied to 
the annual four-cell contingency tables. The procedure tested 
the null hypothesis that the two dichotomous variables are 
independent in the statistical sense, that is, that the proba- 
bility of an observed value of one variable is not affected by 
an observed value of the second variable. Using a 0.95 confidence 
level, the null hypothesis of independence was not rejected 
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for five of the contingency tables: negative lapse rate 
(125 to 80 mb) versus wind shear (350 to 175 mb); negative 
lapse rate (80 to 60 mb) versus wind shear (350 to 175 mb); 
positive lapse rate (175 to 125 mb) versus wind shear 
(700 to 350 mb); positive lapse rate (175 to 125 mb) versus 
wind shear (350 to 175 mb); and positive lapse rate (175 to 125 mb) 
versus windspeed (300 mb). The computed chi-square values for 
most of the other contingency tables were statistically unreason- 
ably high so that the validity of the use of this test for all 
contingency taUt>les was questioned. 

The chi-square test is based on the assvunption that the 
individual observations of a variable are randomly selected 
from the total population of possible independent observations. 

The construction of the present data set negates the assumption 
of independence and thus invalidates the chi-square test. 

Time and space correlations are inherent in the data base, 
since 12-hour sequential observations taken at several stations 
over areas which tend to share homogeneous terrain features 
and synoptic patterns are combined into one set. It is well 
known that dynamic regimes in the stratosphere tend to change 
slowly. Thus, several sequential radiosonde observations over 
a homogeneous horizontal area can actually reflect the same 
atmospheric regime. This correlation between observations 
means that the contingency tables contain multiple counts of 
essentially the same observation. These multiple counts lead 
to the faulty chi-square test conclusion that a high degree of 
dependence exists between variables. A similar problem was 
previously addressed in reference 22, which presents a method 
for estimating the number of independent observations contained 
in the total sample population. Another deficiency of applying 
the chi-square test to annual data samples is the possibility 
that the variables tested may be positively correlated in one 
season but negatively correlated in another. The data set 
used for this pilot study was not designed for such evalua- 
tion; however, the empirical results are relatively informative 
in spite of the lack of rigorously established statistical 
validity. 

Seasonal percent frequencies of the joint occurrence of 
positive and negative lapse rates and high altitude wind 
shears with the other variables meeting the criteria listed in 
table 2 are presented in tables 3, 4, and 5. These tables 
summarize the values obtained from one cell in each contingency 
table that contains the joint occurrence of strong values for 
both variables. Table 3 lists the joint empirical percent 
frequencies of the simultaneous occurrence of positive lapse 
rates equal to or greater than 4** C per kilometer with each of 
the other wind variables listed in table 2. Table 4 is similar 
except that a negative lapse rate (inversion) of less than or 
equal to -5* C per kilometer is the main variable. Table 5 
presents the percent frequencies of the joint occurrence of 
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high wind shears in separate layers and also with high windspeeds 
at several levels. Several other joint combinations of variables 
and layers are possible; the ones depicted herein were subjec- 
tively chosen as being of most interest to high altitude 
turbulence situations affecting supersonic cruise aircraft 
operating over the western United States. 

It should be noted that the use of the term "joint occurrence" 
in this report does not refer to the simultaneous occurrence 
of variables within a given sublayer. The term refers only to 
the simultaneous occurrence of values within ^e main layers. 

For example, a radiosonde observation could yield a maximum 
wind shear between 130 and 120 millibars and a maximum positive 
lapse rate between 105 and 90 millibars. The contingency 
table would contain a count for the joint occurrence of wind 
shear and positive lapse rate within the 125 to 80 millibar 
main layer even though the individual variable values did not 
occur simultaneously within the scune sublayer. 

The joint occurrence of stratospheric strong positive 
lapse rate or reduced buoyant stability and high wind shear is 
associated with low Richardson number and therefore suggests 
an increased probability of turbulence. From table 3 it is 
apparent that this condition is often observed in the summer 
(61 percent of the soundings) in the 175 to 125 millibar layer 
and in more than a third of the soundings in the 125 to 
80 millibar layer. It also occurs almost half the time in 
winter in the 175 to 125 millibar and 125 to 80 millibar 
layers. Note that the wind shear and lapse rate threshold 
values, 0.005 per second and 4® C per kilometer, respectively, 
are equivalent to a Richardson number threshold which is an 
order of magnitude less stringent than the theoretical criteria. 

The joint occurrence of positive lapse rate with high 
values of the wind shear windspeed product does not appear to 
be as strong through 80 millibars. Above 80 millibars, the 
statistics for wind shear and the wind shear windspeed product 
are approximately equivalent. Table 3 also shows that about a 
third of the time positive lapse rates from 175 to 80 millibars 
in winter are coincident with high windspeeds at 700 and 300 
millibars. As altitude increases above 80 millibars in winter, 
the joint occurrence of high lapse rates with the other variables 
decreases. In summer, the joint frequencies are nil. 

Joint percent frequencies for negative lapse rates or 
inversions with other variables are shown in table 4. The 
highest value, 44 percent, occurs in summer for the lapse rate 
wind shear combination in the 125 to 80 millibar layer. This 
combination also occurs in the 80 to 60 millibar layer with 
approximately equal frequency throughout the year. During the 
winter, the joint occurrence of strong inversions and high 
values of each of the other variables, except low altitude 
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wind shears, is observed about 20 percent of the time. In 
summer, most of the variables are jointly observed infrequently. 

TeUole 5 shows the percent frequencies of high altitude 
wind shears occurring with low altitude shears and high windspeeds 
at the selected levels. The frequencies are low in all cases 
in summer. The winter correlations are strongest for wind 
shears in the layers between 175 and 60 millibars and windspeeds 
at 700, 300, and 100 millibars. High 70 millibar windspeeds 
in winter are observed almost half the time that strong shears 
in the 80 to 60 millibar layer occur. At all six levels high 
speeds appear to be observed more often with strong wind 
shears than with either positive lapse rates or inversions. 

Further comparisons were also made between the bivariate 
frequencies observed and those calculated on the assiunption of 
independence. Strong dependence was not exhibited below 
80 millibars. For the layers above 80 millibars, wind shear 
showed some dependence on the windspeed at altitudes of 
100 millibars and higher in all seasons. For excunple, the 
joint occurrence of wind shear in the 60 to 40 millibar layer 
and strong 100 millibar windspeeds in the spring season was 
17 percent, in contrast to less than 11 percent calculated on 
the assumption of independence. Positive lapse rates in the 
layers above 80 millibars indicated some statistical dependence, 
but not as uniformly as wind shear. All joint variables 
showed higher frequencies than calculated (on the assumption 
of independence) by a factor of one- third or more, for at 
least one season in the layers above 80 millibars. Similar 
dependence was not shown for negative lapse rate except in the 
60 to 40 millibar and 40 to 20 millibar layers, where they 
tend to be associated with strong wind shears and windspeeds. 

As an example, the observed joint frequency of negative lapse 
rates and strong wind shears in the 40 to 20 millibar layer 
was approximately twice as great as calculated for all seasons. 

The dependence found between joint variables may be 
sensitive to the nearness of the selected thresholds to physically 
critical values for the layers of interest. 


DISCUSSION 


Application of Climatology to Turbulence Studies 

Climatological information, when available, can benefit 
turbulence forecast studies in three separate areas. First, 
it can allow the investigator to evaluate the meteorological 
parameters which appear to be associated with specific data 
samples. Second, climatology for the appropriate meteorological 
parameters provides a means for estimating seasonal variations 
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in the occurrence of turbulence on the basis of limited amounts 
of experiment:al data. Finally, when predictors are obtained 
from a limited amount of observational data, their performance 
can be evaluated on an annual and seasonal basis by the use of 
climatological information. These applications are discussed 
in the paragraphs below. 

Meteor ological parameter association with turbulence . — ^The 
meteorological conditions observed in association with turbulence 
are usually summarized in terms of the observed frequencies 
within each turbulence intensity category which fulfill specified 
criteria for the magnitude of the meteorological variable. 

The average value of th* meteorological variable within each 
turbulence category is dso used to expose variables which 
influence the turbulence intensity. Comparisons of the results 
from studies of high altitude turbulence (HAT) with statistics 
obtained in the present study are presented below. 

High altitude turbulence results reported in reference 16 
are presented in table 6 in terms of the frequencies of samples 
in each intensity category which exceeded specified meteoro- 
logical criteria. Corresponding frequencies from the present 
sample of upper air data for the western United States are 
also shown on an annual and seasonal basis. For the area of 
^he reference 16 HAT sample, the frequency of maximum wind 
observations exceeding 36 meters per second is believed to be 
closely approximated by the frequency of 300 millibar winds 
exceeding 30 meters per second. Cn this basis, the maximum 
wind values for the turbulence samples exceed the threshold 
criteria for 36 and 37 percent of the cases in the T^ and T^ 

intensity categories and appear to be representative of the annual 
wind background. However, the maximum wind values for the T 2 

intensity samples exceed the threshold criteria in 69 percent 
of the cases; therefore, the frequency for the T 2 category 

differs significantly from the frequency for the other categories 
and the annual background. For this reason, the stronger 
turbulence is believed to be associated with higher speeds in 
the upper troposphere. It is further inferred that it is more 

in the winter season and rare in summer. The frequen— 
the 100 millibar windspeed exceeding 21 meters per second 
exhibit even more contrast. The Tq and category frequencies 

are significantly lower than the annual frequency, but the 
frequency for the T 2 intensity category is significantly greater 

than the winter frequency as well as the annual frequency. 

Unlike the present data, which were separated into 700 to 
350 millibar and 350 to 175 millibar layers, the lower altitude 
wind shear frequencies for the reference 16 HAT data were not 
separated into layers. For comparison the present data are 
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shown for each of these layers, and the frequency for the com- 
bined layers is taken to be less than or equal to their sum. 

The frequencies of strong lower altitude wind shear for the 
Tq and categories are near the winter season frequencies 

indicated by the present data. However, for the T 2 category, 

the frequency of high wind shear is also significantly greater 
than either the annual or seasonal frequencies. 

wind shear, lapse, and inversion rate data for the reference 16 
HAT samples were the maximum values measured within 610 meters of 
the flight altitude in turbulence. Most of these turbulence sauaples 
were obtained from flight altitudes between 12.2 and 21.0 kilometers, 
which may be represented by the present data for layers between 
175 millibars and 40 millibars. These layers for the present data 
are thicker than the ±610 meters for the reference 16 HAT samples, 
and therefore their frequencies are somewhat higher than they would 
have been had the samples been taken from thinner layers. In addition, 
the annual statistics and the seasonal patterns for the maximum wind 
shears, lapse rates, and inversion rates change significantly with 
altitude. Therefore the table contains data for layers which represent 
the range of frequencies and seasonal patterns as well as an overall 
annual average value. 




Wind shears and lapse rates near the HAT altitude equal or 
exceed the specified values at increasingly greater frequencies for 
the Tq, Tj, and T 2 intensity categories. The frequency for strong 

inversion rates is least for the Tj^ category. The frequency for the 


T 2 category is significantly greater than tne background frequency 

indicated by the present data. In consideration of the relatively 
shallow layers used with the HAT data, the frequency of strong 
inversions for the Tq samples is assumed to be at least as great 

as the background frequency. An overall observation of the 
table 6 data is that the specified ranges for all the variables 
are associated with T 2 intensity HAT. Weaker values of the 


variables typify background conditions except for the 100 millibar 

windspeed, Vj^qq, and the flight altitude wind shear, 

which also appear skillful in designating very light turbulence 
or nonturbulent conditions. 


Another study (ref. 19) reported results for similar meteoro- 
logical variables associated with samples obtained by the USAF 
high altitude clear air turbulence (HICAT) program. Data from 
reference 19 and corresponding frequency data from the present 
sample are presented in table 7. The specified ranges used in 
reference 19 are slightly different from the ranges specified 
in reference 16, and a larger altitude interval (±2000 meters 
above and below the flight level) was used to obtain data for 
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the strongest lapse and inversion layers. One effect of this 
procedural difference is a tendency for the frequencies in the 
three HI CAT turbulence intensity categories to be more closely 
grouped near the frequency for all data in the study. Also, 
the frequencies of strong inversions are significantly greater 
in the HICAT sample than in the present data. 

Average values in each intensity category were also deter- 
mined for some of the meteorological variables in the HAT and 
HICAT samples. As shown in tables 8 and 9, the average values 
of the variables generally increase with turbulence intensity. 
Different procedures were used in the two studies, so direct 
comparison of all variables is not feasible. It may be noted 
that the maximum windspeed shows less variation with intensity 
for the HICAT sample. This may be due to the wide variation 
in geographical areas and weather circulation patterns from 
which the HICAT data were obtained. Regional differences in 
the statistics of the meteorological variables may impact the 
experimental results as much as different experimental procedures. 

Both the average values of the meteorological parameters for 
the turbulence intensity categories and the frequencies at hich 
they exceed selected values allow the investigator to interpret 
the circumstances attending the data samples. For example, the 
average value of a variable for the nonturbulent and/or light 
turbulence categories may be near the climatological mean but 
significantly less than the mean for the more intense category. 

In this case it would be assumed that high values of the variable 
would be associated with an increased probability of strong 
turbulence. On the other hand, low values of the variable may 
not appreciably increase the probability of weak or nonturbulent 
conditions. Thus, the availability of the appropriate climatology 
for the meteorological variables would enable investigators to 
judge whether specific variable ranges can discriminate conditions 
for intense turbulence, nonturbulence, or both. 

Turbulence occurrence climatology . — Another application of 
climatological data is the interpetation of results from limited 
turbulence sampling data to more general rates or probabilities 
on a seasonal or annual basis. Such rates can be either turbu- 
lence occurrence rates or turbulence forecast verification rates. 
Both turbulence occurrence rates and the meteorological variables 
in the experimental sample may be influenced by weather pattern 
anomalies existing during the data collection period. To translate 
the results from such data into rates typical of the seasonal or 
annual periods, the data are first examined to determine the 
variables and the magnitudes which most strongly discriminate 
between the turbulence categories of interest. This may be 
accomplished in terms of either conditional probabilities or 
multiple regression techniques. Climatology for the variables 
of interest should then be acquired to determine the probability 
distributions for the magnitude ranges used and to identify 
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' rin^iiv the turbulence occurrence 

dependence between -re^converted to rates for the 

ships or conditional probabilities. 

Estimating forecast verification rates .— Bias 
sa«ples due to size 

-ir^rce” rTndivJusi^^^^^ 

J“«s‘:rbe‘a“oS:p!!^hed"S: llVlUis of thS 

methods of turbulence encounter. P(TR). is used 

for the p^obabil^ that turbulence will be fore- 

n].lt Sr ^aseJ when it is encountered. P(TF/TR). and for cases 
cast loi cases i pitf^nr^ It is assumed that the 

:!:-\!:ertar^r^r2 ueocies foW 

even though P(NR/NF). For convenience of expres- 

aurdrftne'thl'ptoblbuuy that turbulence is not 
encountered as 

P(NR) = 1 * P(TR) 

the probability that turbulence was not forecast even though it 


was encountered as 


P(NF/TR) = 1 - P(TF/TR) 


and the probability that turbulence was not forecast when it was 
not encountered as 

P(NF NR) - I - P(TF NR) 

Kst.^ales for the verification rales desired are then given by: 

P^TF TR) P(TR) 


P(TR/TF) = p(TF/TR) pJtrT ♦ P(TF/NR) P(NR) 

P^NF TR) P(TR) 


P(TR/NF) = P(NF/TR) P(TrT ♦ P(NF/NR) P(NR) 
P(NR TF) - 1 - P(TR/TF) 
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P(NR/NF) = 1 - P(TR/NF) 


Sample Size Considerations 

The primary purpose of the climatology studies discussed 
herein is to establish the statistics for selected meteorological 
variables which are related to the occurrence of clear air 
turbulence. Univariate frequency distributions, conditional 
distributions, bivariate contingencies, conditional probabilities, 
and correlation properties are all applicable to the interpre- 
tation of the meteorological data associated with turbulence. 

The importance of sample size for the determination of the 
distribution functions and statistical confidence levels can 
be appreciated by inspection of the empirical curves and by 
evaluation of statistical criteria. 

Several of the curves obtained for the present sample 
exhibit a lack of stability in their distribution functions. 

This is particularly evident at the tail of the distributions, 
where a limited number of observations are available and 
individual points can greatly modify the shape of the distrib- 
ution curve. The lack of a straight line or a consistent 
simple curvature pattern in the central part of the distribution 
is another indication of insufficient statistical stability. 

Thus, the foremost drawback in this initial study is the small 
size of the sample used. Although the present results only 
give a gross approximation of the statistics, they do provide 
helpful indications of the relative frequencies for the threshold 
magnitudes and seasons of interest. 

The sample size required for the desired degree of statistical 
confidence can be estimated on the assumption that the observa- 
tions are independent and that the distribution is known, 
convenient expressions of statistical confidence are commonly 
available for the mean value and variance of Gaussian distributed 
variables. To illustrate the effect of sample size it will be 
assumed that the measured variables are either Gaussian ot 
have been transformed into variables with a Gaussian distribution. 

Let us specify that it is desired to know the variance, o , to 
within 110 percent and with a certainty of 95 percent. In 

other words, the ratio of the measured variance, s , to the 

true variance, must be between 0.9 and 1.1. To attain 95 

percent confidence that the results are within this interval, 

the probability that the ratio s^/a^ is less than 0.9 wust k>e 
limited to 0.025, and at the same time the probability that the 
ratio is not greater than 1.1 must be at least 0.975. The 
number of independent observations needed is found from the 
inequality 
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m 


1 


^g/2 
n - 1 



^l«g/2 
n - 1 


2 

where X i® the chi-square parameter, a is given by 1 minus 
the confidence level desired, and n is the number of obser- 
vations. The values of the ratios 


and 


X 


0.025 


n - 1 - 


0.9 


are available in reference 23 and indicate sample size require- 
ments of 740 and 800 independent observations, respectively. 
Observation requirements increase rapidly with the accuracy 
desired. For example, 3200 observations would be needed if it 

were desired for the observed variance (s^) to be within ±5 

percent of the true variance (n^) instead of within ±10 percent. 

Sample size requirements for nearness of the sampled mean 
value, X, to the mean value of the population. at specified 

confidence levels are found similarly for Gaussian variables. 

The t-distribut ion is used to relate the desired tolerance on 
the mean value estimate, ^ to the sample variance for 

specified confidence levels where 


xli s 

The desited sample size, n, is found by use of the inequality 

>n - s - \h 

When n exceeds 120, t becomes constant for given levels of 
confidence and the evaluation is analogous to the case of 
sampling a population with known variance where the standardized 
normal distribution variable, z, is used. 


17 


*1-.. * • “x ^ !o/2 

-VS- ' - vn 

For « confidence level tSe“»pU J^'n 

number of ™im“ely 0 1 stendardiied 

value is accurate to «fhrn pp *i observations found for 

unit is 384. If the sample site ot o j,,* ,s percent 

the variance criteria « “«d it would Pr ^jy,in 

certainty that the sample mean value is 
approximately 0.07 standardized unit. 

we can consider Jonf idSncrieJels^'discuslL^*^^*''^ 

required to provide the aiven observation location will yield 

above. If we assume that require at least 10 

an independent ®„^20 stations spaced sufficently to 

years of record ^^?"_tatis?ical inW«dence and stability, 

yield the necessary degiee of aiven as study requirements but 

These sample size technical considerations necessary to 

only to be descriptive of the tecnnicai 

achieve stable results. 

Recommendations 

The present results ^^he^'Lnef itrof''an‘^expanded 

a significantly longer Period of lecor. distributions for 


stud? would be to more f gradie^ wind shear, and so 
the small scale ^®^tical temper gi conditions observed 

forth, and to determine their dependence 

at the standard synoptic levels^ thrLhold, persistence, periodic 
joint occurrence as a -annual variation could be examined 

seasonal oscillation, and study. As noted in the 

but were beyond the scope of the present^^^^ obtained is determined 
previous section, the statistic i degrees of freedom 

by the number of Studies provide results which are 

in the data set. Seveial previous stuai^^^P^^^ 

indicative of temporal and spat i 1 

components. However, ^rvations for other variables, 

degVee of correlation gradient. Therefore, 

such as vertical wind the single station 

the first task ^«^®"®'®"^®^.tions for the variables of interest. 

univariate temporal coiielationsf .delude autocorrelation, 

fp^eitra^'U'Sairov'cbxrtrcrnir- U*f. 24,. 

spatial dependence ^tafer ^^iw'et.“’^b^rrvUronr 

as temporal Spa? ll «lults for autocorrelation and 

substituted for time. ®P^tiai terms of lag distance 

Markov chain techniques wou^^^^ spatial domain, 

wLld^be*??' terms of «avelen,th instead of frequency 
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units (for example, cycles per year), sealo^^ 

and spatial correlation should be accomplished for each season 

as well as on an annual basis. 

Once the univariate correlation characteristics 
and space have been established, an appropriate data set can 
Se created for multivariate analysis Either 
mav be indicated by the results at this point. First, the 
data set may be formed by including only 

are separated sufficiently to insure ade<iuate independence. 

If this is not feasible, more frequent observations may be 
i«d tut tilfnimber of' observations will need to be adjusted 
to reflect the actual number of degrees of freedom (ref. 22) 
for evaluation of statistical confidence in the results. 

Deterministic periodic oscillations which significantly 
rontribute to the total variance may be identified by the 
spectral analyses for some of the variables. ® 

are usually related to the known periodicities of Physical 
phenomena and can be eliminated from the 

yield a new time series comprised mainly of sampling variations 
and noncyclic persistence. 

The thickness of the layers defined in the present study 
was not uniform. Layers were selected for convenience in 
correlating winds reported at the mandatory meteorologica 
analysis JIvels with wind shears and temperature gradients at 
nearbv altitudes. The pressures separating the layers in the 
Diesent study were spaced equally between the mandatory levels. 

?n the sIratLphere the geometrical layer thickness ranged 

frf.m Tflirmete^rfor the layer between 80 millibars and 60 millibars 

to 4481 meters for the layer between 40 emoiriral^^^ 

Since a uniform layer thickness was not used, the empirical 
urSblbtliTies foi Lch layer are biased by thickness. Therefore 
it is recommended that a uniform layer thickness be defined and 
incorporated in future studies. 


CONCLUDING REMARKS 


A Pilot study was performed to characterize the statistical 
behavior of the maximum vertical wind shears and 
gradients occurring in specified altitude layeis^ The y 

ii«ied data from 20 upper an stations in the western Unite 
s* ^tes tor a i year period of record and emphasized layers in 
thf lowe? sUatospU?e. In general, strong wind shears, lapse 
atef Ind inCe?lions wee observed less frequently as altitude 
inctlaseS frl l^S millibars to 20 millibars. 
basis, the shears and gradients were *" "‘J'" *2 

weaker in summer, but minor deviations to the P^5^^®^” ^®^® 
no?e5 in association with increased tropopause altitude in 
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summer and with the stratospheric wind reversal in the spring 
and fall. 


The empirical probability distribution curves tended to 
follow a log-normal distribution curve. Because of the limited 
sample size, the distribution curves showed indications of in- 
stability, but even so, they provide a coarse estimation of 
the percentile frequencies for the higher magnitudes of interest. 

The joint occurrence of variables exceeding specified 
magnitudes was ex6unined with bivariate contingency t^les. 

The contrast between univariate probability and conditional 
probability for several pairs of variables indicated weak 
dependence. However, these cases of dependence were essentially 
limited to altitudes and seasons when the frequencies of the 
strong magnitudes were relatively low. For most of the contin- 
gency tables the computed chi-square values were unreasonably 
high. This is believed to be due to serial correlation between 
samples for each of the variables, and therefore valid conclusions 
on the degree of dependence between variables could not be 
established with the present chi-square test results. 

Applications of the univariate and bivariate results to 
studies of turbulence encountered by aircraft were discussed. 
Examples presented included the translation of limited aircraft 
turbulence encounter experience into the amount of turbulence 
and the forecast verification rates which could be expected on 
a climatological basis. 

Recommended steps for improving the present results are 
given below. 

1. Extend the data base to a significantly longer period of 
record. 

2. Analyze univariate temporal correlation (dependence) for 
geographically representative single station records. 

3. Analyze univariate spatial correlation. 

4. Evaluate the period of record, number of stations, and the 
spacing of observations needed to provide the desired statistical 
confidence . 

5. Establish appropriate data sets and implement multiple station 
multivariate analyses. 

Dryden Flight Research Center 

National Aeronautics and Space Administration 
Edwards, Calif. 93523, Nov. 12, 1900 
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CLASS INTERVALS, AND RESTRICTIONS 















































SEASONAL EMPIRICAL BIVARIATE PERCENT FREQUENCIE: 
■OR VARIABLES MEETING CRITERIA IN TABLE 2 
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TABLE 8.— AVERAGE VALUES OF METEOROLOGICAL VARIABLES 
FOR THE RJSFEREKCE 18 TURBULENCE INTENSITY CATEGORIES 


Variable 


All data 



m/ufsc 

34 

'^lOO' 

m/sec 

12 


per sec 

.012 


TABLE 9. — AVERAGE VALUES OF METEOROLOGICAL VARIABLES 
FOR THE REFERENCE 1> TURBULENCE INTENSITY CATEGORIES 


Variable 




•C/kiR 


^ Li ght 


Light 


^.Light 
to moderate 

All data 

42 

38 

S.2 

3.4 
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6.3 
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Flvlure winter cumulative distributions of wind speeds 
for Indicated pressure levels. 
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Fidure 3. Spring cumulative distributions of wind speeds 
for Indicated pressure levels. 
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Figure 6. Seasonal cumulative distributions of 700 mb wind speeds. 
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Figure 9. Seasonal cumulative distributions of 70 mb wind speeds. 



Wind speed, m/sec 


Figure 10. Seasonal cumulative distributions of 50 mb wind speeds. 


















wind shear. 1/soc 


Moure 18. Seasonal cumulative distributions of wind shears 
^ for the 700-350 mb layer. 
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Figure 23. Sprin^j ciimuKitivo Jlstr ihut Ions of VNin.l sne.i»*s 

for i»Hlkrttod layers. 
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Fiijure 24, Su'vier cj'’Uilativc J i stributlor.s of wind shiMrs 
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Figure 25. Fall cumulative distributions of wind shears 

for indicated layers. 
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for 'nlicated layers. 


71 






r 

E 

I 

I 


fe 

I 


f 

i 

I 




Frequency 

percent 





Pressure, 

mb 

~ 40-2Q 

- - 90-m 

80-60 

t28-80 

175-125 

350-175 

700-350 

— — Troiwpeuse 
to 20 


c i.i L- 
6 7 8 9X 10’ 




^'ositive Itipsr r*Ue. 

Miuf t' 4i. .iinfor L cir\j],it ive J i s i r ihut Ions of positive lapse rates 

for induated layers. 


72 


ORIGINAL PAGE IB 
OF P«X)R QUALIIY 
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Fiqure 44, Spring cumulative distributions of positive lapse rates 

for indicated layers. 
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Pi'situt' lapse rate, C/kn! 


fiqurc 4S. Sa-ner arn.lative disi ri^Jtions of positive lapse rates 

for Indicated layers. 
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Fiaure 46. fall cu’ulative ilisl'ibutions of positive lapse rates 

for indicated layers. 
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Posit We lapse rate, *C/liin 


Fiaure 47. Seasonal cumulative distributions of positive lapse rates 
^ for the 700-35C mb layer. 











figure 51. Seasonal 


cumulative distributions of positive lapse rates 
for the 30-60 mb layer. 
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M:iuro i)3. Se<^sonal Oi-ulative distributions of positive lapse rates 

for the 40-20 nb layer. 
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F i ijre »»inter cunuldtive distributions of inversion rates 
(ntn^ative lapse rates) for Indicated layers. 
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Inversion rate, "C/Km 


'■Kjure S? . iprinj cunulattve distributions of Inversion rates 
(ne>^at1ve lapse rates' for Indicated layers. 
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59. F.^’l .-..•n.litive distributions of Inversion rates 

rdt^s) for ind^cdt.od layers. 













